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PREFACE 

This  paper  is  one  of  a series  of  IDA  papers  on  combat 
modeling,  sponsored  by  the  IDA  Independent  Research  Program. 
Mathematical  attrition  processes  are  fundamental  to  the  com- 
bat simulation  models  used  In  many  studies  of  defense  problems. 
The  DOAE  research  on  attrition  processes  is  a significant  body 
of  work,  and  the  comments  of  this  paper  are  Intended  to  present 
and  evaluate  the  DOAE  results  in  the  context  of  other  mathe- 
matical, computational  and  applied  work  in  combat  modeling. 
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1.  INTRODUCTION 

This  paper  is  a review  and  summary  of  the  following  memo- 
randa prepared  by  the  Defence  Operational  Analysis  Establishment 
(DOAE)  of  the  UK: 

(1)  The  Mathematics  of  Battle  I:  A Bivariate  Probability 
Distribution,  by  T.G.  Weale; 

(2)  The  Mathematics  of  Battle  II:  The  Moments  of  the 
Distribution  of  Battle  States,  by  T.G.  Weale; 

(3)  The  Mathematics  of  Battle  III:  Approximate  Moments 
of  the  Distribution  of  States  of  a Simple  Hetero- 
geneous Battle,  by  N.  Jennings; 

(4)  The  Mathematics  of  Battle  IV:  Stochastic  "Linear 
Law"  Battles,  by  N.  Jennings; 

(5)  The  Mathematics  of  Battle-V:-  Homogeneous  Battles 
with  General  Attrition  Functions,  by  T.G.  Weale; 

(6)  The  Mathematics  of  Battle  VI:  The  Distribution  of 
the  Duration  of  Battle,  by  T.G.  Weale; 

(7)  The  Mathematics  of  Battle  VII:  Moments  of  the  Dis- 
tribution of  States  for  a Battle  with  General 
Attrition  Functions,  by  T.G.  Weale  and  E.  Peryer. 

These  papers  are  hereafter  cited  as  references  [13],  [1^1,  [5], 
[6],  [15],  [l6],  and  [17],  respectively. 

Common  to  all  seven  papers  is  an  attempt  to  deal  on  a 
computational  basis  with  certain  stochastic  attrition  processes 
of  the  sort  discussed  in  [10].  Not  all  processes  discussed  in 
the  DOAE  memoranda  are  mentioned  explicitly  in  [10];  Indeed  some 
of  the  work  in  [15,  l6,  17],  wherein  essentially  arbitrary 
attrition  functions  are  oermltted,  is  not  physically  justified. 
There  is  no  doubt  that  Markov  processes  exist  with  the  genera- 
tors Indicated  (cf.  Section  6 for  details);  whether  such 
processes  correspond  to  physically  plausible  (or  even 
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physically  definable)  sets  of  assumptions  Is  uncertain.  Further 
remarks  concerning  this  point,  which  Is  of  considerable  Impor- 
tance, may  be  found  In  Sections  6 and  9;  cf.  also  [10]. 

The  emphasis  In  this  review  Is  on  assumptions  underlying 
the  attrition  processes  discussed,  on  mathematical  computations 
and  approximations  presented  In  the  papers,  and  on  probabilis- 
tic derivations  and  Interpretations  of  certain  results  that  have 
been  derived  In  the  DOAE  papers  using  methods  from  the  field  of 
different  equations.  In  particular,  we  have  not  reviewed  care- 
fully the  computer  programs  Included  In  the  papers,  nor  have  we 
analyzed  either  sample  outputs  Included  In  the  papers  or  (except 
In  preparation  of  the  related  papers  [8]  and  [9])  a large  number 
of  additional  outputs  generously  provided  to  us  by  DOAE.  The 
few  analyses  performed  have  shown  DOAE  results  to  be  consistent 
with  those  the  present  author  reported  In  [8]  and  [9]* 

We  are  aware  that  the  main  purpose  of  the  DOAE  work  was  to 
produce  computer  programs  with  which  numerical  conclusions  could 
be  obtained,  and  that  much  of  our  criticism  Is  directed  at  other 
aspects  of  the  work.  Even  so.  It  Is  useful  to  have  available  a 
careful  discussion  of  the  foundations  of  the  processes  under 
analysis,  which  Is  what  we  have  attempted  to  provide. 

We  wish,  at  the  beginning  of  this  review,  to  commend  the 
authors  of  all  seven  DOAE  papers  for  uniformly  high  technical 
and  expository  qualities  of  their  work.  The  papers  were  a 
pleasure  to  read. 

Throughout  this  paper  our  notation  and  terminology 
concerning  Markov  attrition  processes  are  those  of  [10]. 


2.  REVIEW  OF  "A  BIVARIATE  PROBABILITY  DISTRIBUTION 


This  paper  reports  efforts  to  compute  and  approximate  the 
transition  function  of  a stochastic  attrition  process  that  is 
analogous  to  the  homogeneous  Lanchester  square  differential 
equations  of  combat.  That  is,  the  author  deals  with  Homogene- 
ous Process  1 of  [lu,  p.  17],  which  has  independent  engagement 
initiation  and  single  kills.  Referring  to  [10],  we  note  that 
this  Markov  attrition  process,  which  is  denoted  by  ( ,R^ ) , 
has  infinitesimal  generator  A given  by  ” 

A((l,J),(l,j-l))  = icg 

(1)  A((i,j),(i,J))  = - (iCg+JCp) 

A((i,J),(i-l,j))  = jCj^  , 

where  Co,c„  are  positive  constants,  jump  function  X given  by 

D n 

(2)  X(l,j)  = ICg  + jCg  , 

and  transition  matrix  Q of  the  embedded  Markov  chain  given  by 


(3) 


ic 


B 


Q((1,J),(1.J-1))  = 15— 


Jc 


R 


S((1,J).(1-1,J))  - tcg  Jc„ 


Rather  than  allow  the  combat  to  continue  until  one  side 
is  annihilated,  the  author  imposes  termination  levels  m^ 
for  Blue  and  mj,  for  Red.  That  is,  all  states  of  the  forms 

A 

(i,mg)  and  (mg,J)  are  absorbing.  For  such  states  a, expressions 
(1),  (2),  and  (3)  are  not  valid  and  one  has.  Instead,  that  for 
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each  state  g 


A(a,B)  = 0 
X(a)  = 0 

Q(a,g)  = I(a,8), 


where  I is  the  Identity  matrix.  Observe  that  the  state  (m„,m„) 

D K 

is  almost  surely  never  entered. 


Let  us  denote  by  (P^)  the  transition  function  of  this 
attrition  process;  cf.  [10,  p.lO].  In  [13] j Weale  is  concerned 
mainly  with  computation  and  approximation  of  the  transition 
function  (P^)j  that  is,  computation  and  approximation,  for  each 
fixed  t,  of  the  Markov  matrix  . Another  interest  is  computa- 
tion of  the  limit  P = lim  P , which  obviously  exists  and  is 

t->-«>  ^ 

the  distribution  of  the  terminal  state  of  the  process. 


Weale  begins  by  deriving,  in  a heuristic  manner  based  on  a 
physical  interpretation  of  the  attrition  coefficients  c_  and 

D 

Cpj,  the  forward  equation  for  the  process,  namely 
(i|)  P^((l,J),(k,)l))  = Cj^JlP^((i,J),(k+l,<l)) 

+ CgkP^((l,J),(k,Jl+l)) 

- (Cgk+Cg«,)P^((i,J)  ,(k,t)) 

= P^A(  (l,j  ) , (k,ll) ) . 

The  expression  (4)  is  valid  if  (k,s,)  is  not  an  absorbing  state; 
for  absorbing  states  (k,mg)  we  have 

(5a)  P^( (i, j ) , (k,mg) ) = CgkP^( (i, j ) , (k,mg+l) ) , 

while  for  states  of  the  form  (mg,*.) 

(5b)  P'((l,J),(mg,)l))  = Cg)lP^((l,J),(mg+l,£))  . 

That  these  equations  are  valid  is  an  immediate  and  rigorous 
consequence  of  the  probabilistic  derivation  appearing  in  [10]; 
cf.  also  the  Appendix  to  [7],  where  a detailed  proof  is  given. 
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The  author  proceeds  to  note  that  for  fixed  Initial  condi- 
tions (i,J)  there  apparently  exists  no  simple  closed-form  solu- 
tion to  the  forward  equations  (4)-(5).  To  the  reviewer's 
knowledge  this  is  still  so,  although  we  note  below  some  steps 
in  this  direction  hinted  at,  but  not  fully  Identified,  by  Weale. 
That  an  explicit  solution  to  (4)-(5)  is  hard  to  obtain  is  not, 
of  course,  a new  piece  of  knowledge.  R.N.  Snow,  who  originally 
proposed  this  particular  model  [12],  was  aware  of  the  difficulty. 

With  (i,j)  fixed,  certain  of  the  forward  equations  can  be 
solved  in  closed  form.  For  example,  it  is  Immediate  that 


P^((i,J)>(iJ))  = e 


-(Cgi+CpjJ  )t 


this  follows  from  (2)  and  properties  of  continuous  time  Markov 
processes  with  finite  state  space;  cf.  [1]  and  [10,  p.ll].  The 
author  gives  a nonprobablllstic  derivation  from  the  forward 
equation,  using  differential  equation  methods. 


Also,  the  author  obtains  the  following  explicit  solutions 
by  nonprobablllstic  reasoning.  For  1 J > 


(6a)  P^((i, ]),(!,£))  = 


U 


1 


-(Cgi+Cp«,)t  / -c^t'J”^ 


while  for  m^  < k 5.  i» 


(6b)  P^((l,j),(k,J))  = ( 


i-k)!\cB; 


-(Cgk+CpjJ  )t 


C. \i-k 


These  equations,  which  correspond  to  one  side's  having  suffered 
no  casualties,  do  not  seem  to  have  appeared  elsewhere.  Weale 's 
derivation,  as  noted  above,  is  based  on  properties  of  the  forward 
equations  as  differential  equations,  so  it  seems  worthwhile  to 
sketch  a probabilistic  derivation.  For  fixed  initial  conditions 
(1,J)  and  a fixed  time  t 
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(7)  P^((l,J),(l,Jl)) 


_ ^ ^ j first  J-H  transitions  are  kills  of 

Reds,  transition  occurs 

s t 

before  t,  (J-Jl+1)  transition 
occurs  after  t} 

= P ^ ^ ^ I first  J-H  transitions  are  kills 
of  Red} 

X jiime  for  (J-A)  transitions  is 

£ t,  time  for  (J-)l+l)  transi- 
tions is  > 1 1 first  (J-H)  transi- 
tions are  kills  of  Reds} 


P=il+1 


■f  <V*=rJ> 


-(Cgi+CRj )t^ 


^t-ti 

j (Cgi+Cj^(J-l)  )e 


-(Cgl+Cp( j-1) )t2 


ip  -(c  l+c_(J-2))t_ 

(C3i+cp(j-2))e  ^ ^ 


.t-Uj-A  + 1 -(c  i + c (ll+l))t 

) c 0+1  '1  'lo  “ ^ 


( Cgi+CpC  Jl+1)  )e 

-(Cgi+CpJl)  (t-Uj_j^) 


where  is  the  probability  law  of  the  attrition  process  subject 

to  ^{(BQ,RQ)  = (l,j)}  = 1 and  Up  = (t^  +...+  tp). 
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In  arriving  at  this  expression  the  first  terms  are  obtained 
from  the  transition  matrix  Q of  the  embedded  Markov  chain  and 
the  integral  term  arises  from  the  Jump  function  X.  One  must 
invoke  the  characterization  given  in  Corollary  (8.3.11)  of  [1]. 
Tedious  but  straightforward  calculations  transform  (7)  to  (6). 

The  author  takes  note  of  the  further  facts  that  for  (k,i!,) 
not  an  absorbing  state  there  exists  the  representation 


1 j -(c  p+Cpa)t 

(8)  P.((i,J),(k,Jl))  = I I c((i,j),(k,il),(p,a))e  ® ” 

p=k  a=Jl 

where  the  c ( (i, J ) , (k, A) , ( p ,o) ) are  constants,  and  that,  for 

mg  < k £ 1, 


(9)  P^((i,J),(k,mg))  = P^((i,J),(k,mg))  + 


I I c(  (i,J  ) , (k,mg),  (p,a))n 


-(CgP+Cga)t 


p=k  a=mg+l 


with  an  analogous  expression  holding  for  states  of  the  form  (mg,*.). 
These  expressions  are  derived  from  general  considerations  about 
solutions  of  systems  of  differential  equations.  If  the  constants 
appearing  in  these  expressions  were  known  in  closed  form,  compu- 
tlonal  problems  would  be  solved.  Weale  asserts  [13,  p.ll]  that 
the  constants  can  be  calculated  recursively,  but  does  not  elabo- 
rate; we  sketch  below  a possible  way  of  performing  this  computation. 
For  ease  of  exposition  fix  (1,J)  and  consider  (k,*,)  = (i,j-l). 

We  need  to  calculate  c ( ( i , J ) , (i , j-1 ) , ( i , j ) ) and  c( (1,J ) , (i,J-l) , 
(i,j-l)),  assuming  that  (i,j-l)  is  not  an  absorbing  state.  Dif- 
ferentiating (8)  and  setting  t = 0 yields,  on  the  basis  of  the 
fundamental  relation  A = P^,  the  equation 


(10a) 


A((i,j),(l,J-l))=  -(Cgi+Cgj)c((i,J),(i,J-l),(i,j)) 


-(Cgl+Cg(j-l))c((l,J),(i,j-l),(l,J-l))  , 


while  setting  t = 0 in  (8)  gives 

(10b)  c((l,J),(i,j-l),(l,j))  + c((l,J),(i,J-l),(i,J-l)) 

= PQ((i,J),(l,j-l)) 

= 0 . 

One  may  then  solve  (10)  for  c ( ( 1 , J ) , ( i , J -1 ) , ( 1 , J ) ) and 
c( (IjJ ) > (1> j“l) » (1>J-1) ) • More  equations  are  necessary  when 
(l,j)  and  (kjA)  are  not  adjacent  states;  exactly  how  these  are 
to  be  obtained  is  not  clear  and  is  a problem  worthy  of  further 
research. 

Concerning  computation  of  the  limit  (and  terminal)  distri- 
bution P the  author  obtains  the  recursion  relations 
”, 

c k 

(lla)  P^(  (i,  j.) ; (k,mp) ) - ^ ^ + c (m^j+l)  ^ ^ ^ » J ^ 5 ( ) ) 

B R R 

(lib)  P.((l,J);(me.»))  = 5-ti,  • 

B B R 

which  are  valid  for  m^  < k < and  m^  < < j , where  the  asterisk 

in  the  right-hand  side  of  (11)  Implies  that  the  probabilities  in 
question  are  computed  with  respect  to  termination  levels  (mcm^+l) 
in  (lla)  and  (mg+l,mj^)  in  (11b),  respectively.  The  result  in 
(11)  is  not  a closed-form  solution  but  is  certainly  feasible  for 
numerical  computations  even  though  to  compute  the  termination 
distribution  for  one  set  of  termination  levels,  the  termination 
distribution  must  first  be  computed  for  every  set  of  higher 
termination  levels. 

Two  alternatives  are  available  for  performing  numerical 
computations.  First,  one  can  introduce  transforms  in  the 
following  manner;  define  the  resolvent  (8^)j^>q  of  the  transi- 
tion function  (P^)  by 


lawiia^^TuartaaaiMMiiiiii 


where  matrix  integration  is  performed  componentwise.  The  matrix 
is  called  the  X-potential  matrix  of  the  transition  function 
(P  ).  Then  either  directly  from  (8)-(9)  or  by  virtue  of  a 
standard  Tauberlan  theorem  [4,  p.  421]  it  follows  that 

(12)  P = lim  XR,  . 

“ X|0  ^ 

It  is  known  [1,2]  that  for  each  X > 0 the  X-potential  matrix  R^^ 

is  given  explicitly  by 

(13)  = (XI-A)"^ 

and  the  Infinitesimal  generator  A is  known  in  closed  form,  so  by 

(12)  for  small  X,  XR.  is  a good  approximation  to  P . Second, 

A 

since  all  recurrent  states  are  absorbing, 

(14)  P^  = lim  , 

n 

where  Q is  the  transition  matrix  of  the  embedded  Markov  chain. 
Hence  large  powers  of  Q are  good  approximations  to  and. 
Indeed,  the  convergence  in  (14)  takes  place  geometrically  fast. 
One  could  thus  compute,  for  example,  Q^,Q^  = (Q^xO.^),Q 
to  quickly  obtain  an  approximation  to  q”  = P^.  Here,  moreover, 
the  computational  work  is  less  than  for  the  method  based  on 

(13) ,  for  the  latter  requires  a time-consuming  matrix  inversion, 
while  the  method  based  on  (14)  requires  relatively  few  matrix 
multiplications.  Neither  of  the  above  methods  is,  however, 
numerically  exact,  whereas  the  method  using  (11)  is. 

The  paper  also  deals  with  numerical  solution  of  the  forward 
equations  (4)-(5);  this  is  a further  useful  contribution  and  we 
believe  that  the  associated  computer  program  has  received  less 
attention  than  it  deserves.  However,  the  usual  difficulties 
of  extrapolating  numerical  results  and  of  extracting  theoretical 


insight  from  empirical  data  remain.  Nonetheless  these  programs 
seem  to  represent  a significant  contribution  to  qualitative 
understanding  of  this  stochastic  attrition  process. 


=1 
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3.  REVIEW  OF  "THE  MOMENTS  OF  THE  DISTRIBUTION 
OF  BATTLE  STATES" 

The  principal  objective  of  this  paper  [14]  by  T.G.  Weale 
is  to  study  time-dependent  behavior  of  moments  of  the  stochas- 
tic attrition  process  ((B^,R^))  with  generator  A given  by 

A((l,J),(i,J-l))  = Cgi 

(15)  A((1,J),(1,J))  = - (Cgi+CgJ) 

A((i,J),(l-l,j))  = CgJ  ; 

' i.e.,  the  same  homogeneous  process  discussed  in  [13] > which  was 

reviewed  in  the  preceding  section. 

2 

For  initial  states  (l,j),  points  (x,y)  f R , nonnegative 
Integers  r,s,  and  times  t ^ 0,  let 

(16)  M^^j(r,s(x,y);t)  = ^[(B^-x)'’(R^-y)®] 

i J r.  « 

j = I I P ((i,J),(k,«,))(k-x)  («,-y)®  , 

k=mg  K'=nig 

where  (Po.)  is  the  transition  function  of  the  attrition  process 
and  m„,mo  are  the  termination  thresholds  for  Blue  and  Red, 
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respectively.  If  the  transition  function  (P^)  were  known 
explicitly  or  a good  numerical  approximation  were  available, 
then  (16)  could  be  used  directly  to  compute  all  moments  desired. 
In  the  absence  of  a closed-form  expression  for  the  transition 
function,  it  is  useful — from  both  theoretical  and  practical 
points  of  view — to  approximate  moments,  to  derive  relations 
among  moments,  and  to  seek  qualitative  information  about 
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moments  by  studying  (16)  directly.  Such  is  the  author's  main 
goal  and  contribution;  we  now  proceed  to  describe  his  work  in 
more  detail. 

Assuming  that  (x,y)  is  not  a function  of  t,  one  may 
differentiate  (16)  to  obtain 

(17)  M;  (r,s,(x,y);t)  = I I P)  ( (i,J  ) , (k,£) ) (k-x)^(il-y 

k A ^ 

= I I P.A((l,J),(k,A))(k-x)^(A-y)® 
k A ^ 

= I I AP. ((i,J),(k,A))(k-x)^(A-y)®  , 
k A ^ 


where  A is  given  by  (15)  and  where  the  second  and  third  equali- 
ties in  (17)  hold  by  virtue  of  the  forward  equation  (^)  and  the 
related  backward  equation  (cf.[l]),  respectively.  Direct  sub- 
stitution for  A in  either  of  the  latter  two  equalities  in  (17) 
does  not  lead  to  a tractable  system  of  equations.  Indeed,  use 
of  AP.  produces  a system  that  is  not  closed  with  respect  to  the 
initial  conditions  (i,J);  l.e.,  solutions  for  (i,J+l)  and 
(i+l,J)  are  needed  in  order  to  obtain  the  solution  for  (1,J). 
Use  of  P^A  leads  to  a closed  system  that,  however,  appears 
very  difficult  to  solve;  in  particular  one  does  not  obtain 
a differential  equation  for  M.  ,(r,s,  (x,y) ; • ) . The  author  clr- 

1 5 J 

cumvents  these  obstacles  by  introducing  the  truncated  moments 

(18)  j(r,s,(x,y);t)  = ^ [ (B^-x)^(R^-y)® ; lT>t } ] 

where  T is  the  random  time  of  termination  of  the  battle,  at 
which  either  = m^  and  > mj^  or  B^  > m^  and  R^  = m^^. 

The  author  obtains  the  following  equation  (equation  (12)  of 
[14]),  which  is  his  main  theoretical  result: 
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(19)  j(r,s,  (x,y)  jt) 


- c. 


I (-l)^”^(n)CW  (r-p,s+l,(x,y) ;t) + 

p=l  P 


yMijvi’-p,s,(x,y)  ;t)] 


If  r = 0 or  s = 0 the  corresponding  summation  (taken  over  an  empty 
index  set)  is  zero.  Except  as  a means  of  deriving  relations  among 
moments  and  truncated  moments,  (19)  seems  to  have  little  applica- 
tion. From  the  standpoint  of  numerical  calculations,  truncated 
moments  are  rather  difficult  to  obtain  and  even  then  (19)  involves 
truncated  moments  of  all  lower  orders  and  produces  only  the  deriva- 
tive M* . It  might  be  better — if  one's  objective  were  numerical 
results — first  to  produce  numerical  approximations  to  the  transi- 
tion function  (P^)  and  then  to  employ  (16)  directly. 

This  leaves  the  possibility  that  special  cases  of  (19)  may  lead 
to  relations  among  the  moments  M and  truncated  moments  M that  extend 
one's  understanding  of  the  attrition  process — at  least  in  a qualita- 
tive sense.  The  author,  therefore,  proceeds  to  consider  several 
special  cases. 

Taking  r=l,  s=0,  x=0,  y=0  leads  to  the  equation 

(20)  = - CpE^^’J^[R^;{T>t}]  . 

This  famous  equation,  derived  originally  by  R.N.  Snow  [12]  states 
that  for  small  values  of  t the  expectations  of  the  stochastic  attri- 
tion process  approximately  satisfy  Lanchester's  square-law  differ- 
ential model  of  combat,  since  an  analogous  expression  is  valid  for 
E(i>J)[R^], 

To  our  knowledge,  a completely  probabilistic  derivation  of  (20) 
has  not  been  given  before.  The  Theorem  below  is  of  interest,  there- 
fore, not  only  for  its  novelty,  but  also  because  it  both  simplifies 
and  extends  the  results  obtained  by  Weale  in  [l4]. 
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(21)  THEOREM.  Let  f be  a function  on  the  state  space  E of  the 
attrition  process.  Then  for  each  t 


(22)  ^E^^’J^[f(B^,R^)]  = E^^’J^[Af(B^,R^);  jT>t|]  , 


where  A is  the  infinitesimal  generator  given  by  (15). 
PROOF.  By  virtue  of  the  forward  equation  (4) 


(23)  ^ E^^’'^^[f(B^,R^)]  = ^ I P^((l,J),a)f(a) 


= I (|t 


= I P. A((i,j);a)f(a) 


= I I P. ((l,J);6)A(B,a)f(a) 
a 6 ^ 


= I P.((i,J);B)  I A(6,a)f(a) 
6 ^ 


a 


= I P. ((i,J);6)Af(3) 
B ^ 


= ^[Af(B^,R^)]  . 


But  since  Af(B)  = 0 when  B is  an  absorbing  state, 

(24)  ^[Af(B^,R^)]  = E^^’‘^^[Af(B^,R^);{T>t}]  . 

The  proof  follows  by  combination  of  (23)  and  (24).  |] 


The  reason  for  using  the  right-hand  side  of  (24)  in  (22), 
rather  than  the  left-hand  side,  is  that  in  particular  cases  the 
function  Af  will  have  a closed-form  expression  for  nonabsorb- 
ing states  which  is  Invalid  for  absorbing  states.  Use  of  the 
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expression  in  the  right-hand  side  of  (24)  is  permissible,  but 
one  cannot  do  this  on  the  left-hand  side  of  (24).  The  examples 
below  illustrate. 


EXAMPLE.  For  f(l,J)  = 1 one  has  for  nonabsorbing  (1,J), 

Af(i,J)  = - 

and  (22)  becomes  (20).  Similarly,  for  f(l,J)  = J,  one  has 


Af(l,J)  = - icg 


and  obtains  the  corresponding  equation  for  ^ 


EXAMPLE.  For  f(i,J)  = IJ  we  have,  provided  (1,J)  is  not 
absorbing, 

Af(l,J)  = - Cgi^  - Cj^J^ 


and  (22)  becomes 

d_  j T>t | ] 

- Cj^E^^”5^[R^;  |T>th  . 

EXAMPLE.  For  f(i,j)  = i^  and  (i,J)  nonabsorbingj 
Af(l,j)  = - 2Cj^lj  + 


and  one  obtains  the  equation 

d_  = - 2Cj^  E^^’J^CB^R^;  lT>tp 

+ jT>tf]  . 

One  may  then  derive  expressions  for  derivatives  of  variances, 
and  so  on. 
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We  remark — and  this  observation  Is  of  some  Importance — 
that  (22)  Is  valid  If  T Is  replaced  by  any  stopping  time  S such 
that  (B^,R^)  = (BgjRg)  on  {S  < t}.  S need  not  be  defined  In 
terms  of  absolute  thresholds,  for  example. 

One  can  also  use  Theorem  (21)  to  derive  the  general  equa- 
tion (19)  as  well  as  expressions  for  derivatives  of  functions 

I*  s 

not  expressible  In  the  form  f(l,J)  = (1-x)  (J-y)  . We  plan  to 
report  more  fully  on  consequences  of  our  Theorem  In  a forth- 
coming paper  [11]. 

Most  of  the  remainder  of  the  memorandum  [1^]  Is  occupied 
with  computations  similar  to  (19).  Many  of  the  author's  results 
are  worked  out  In  painstaking  detail  and  In  this  respect  he  has 
made  a meaningful  contribution.  We  know  of  no  other  reference 
for  many  of  these  calculations. 

In  the  very  last  part  of  the  paper,  Weale  performs  a simi- 
lar analysis  of  an  analogous  heterogeneous  attrition  process. 

The  process  considered  Is  similar  but  not  Identical  to  Hetero- 
geneous Process  1 of  [10],  with  Independent  engagement  Initiation 
and  single  kills.  If  there  are  M weapon  types  on  the  Blue  side 
and  N weapon  types  on  the  Red  side  and  If  the  states  of  the 
process  are  represented  as  vectors  of  the  form 

(x,y)  = (x^, . . . ,Xn^;  y^, . . . ,yjj)  , 

then  the  Infinitesimal  generator  A of  the  associated  (vector- 
valued) attrition  process  ((B.  ,R.  ))  Is  given  by 

U 1/ 

M 

A ( ( X ,y ) , (x;  y^ , . . . ,yj -1 , . . . ,yjg) ) - Cg(l,j)x^ 

M N 

(25)  A(  (x,y) , (x,y) ) = - I I (og(l,j)x  + Cg(.1,l)yj) 

1=1  J=1  in  J 

N 

A(  (x,y) , (x^, . . . ,x^-l, . . . ,Xj^;y) ) = ‘ 
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The  reader  is  referred  to  [7]  and  [10]  for  details  concerning 
this  process.  Termination  states  are,  not  entirely  reasonably, 
taken  to  be  states  (x,y)  for  which  either 


(26a) 


M 


1=1 


“1^1  1 '"b 


or 


(26b) 


- "r 


where  the  a.,  the  3,,  m_  and  are  prescribed  in  advance.  In 
1 j B K 

such  states  boundary  conditions  of  the  form  (5)  are  valid,  rather 
than  (25).  Certain  modifications  are  necessary  because  some  of 
the  or  y^  may  be  zero  without  either  (26a)  or  (26b)  being  satis- 
fied. The  author  then  notes  the  forward  equation 


for  the  transition  function  (P^)  of  this  attrition  process, 
states  (seemingly  correctly)  that  it  cannot  easily  be  solved  in 
closed  form,  and  remarks  about  its  general  form;  cf.  equations 
(8)  and  (9)  and  associated  comments  in  Section  2. 

Finally,  he  derives  equations  analogous  to  (19)  for  various 
moments.  Appropriately  extended.  Theorem  (21)  may  more  easily 
yield  the  same  results.  Nonetheless,  the  equations  presented 
stand  as  a contribution  because  of  their  explicitness  and  the 
care  used  in  their  derivation. 
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REVIEW  OF  "APPROXIMATE  MOMENTS  OF  THE  DISTRIBUTION 
OF  STATES  OF  A SIMPLE  HETEROGENEOUS  BATTLE" 


This  paper  [5]  by  N.  Jennings  is  not  of  the  same  Interest 
or  Importance  as  the  papers  reviewed  in  the  two  preceding  sec- 
tions; much  of  it,  indeed,  concerns  calculation  of  error  bounds 
whose  theoretical  and  practical  value  is  uncertain.  In  addition 
the  author  seems  not  to  distinguish  relatively  Important  ideas 
from  essentially  uninteresting  (albeit  involved  and  complicated) 
computations . 


Consider  the  heterogeneous  stochastic  attrition  process 
with  the  generator  A given  by  expression  (25)  in  Section  3; 
this  process  is  discussed  in  more  detail  there  and  in  [7,10]. 
Instead  of  the  unrealistic  termination  rule  embodied  in  (26), 
the  author  chooses  the  following  rule:  for  each  1 = 1,  ...,  ! 
there  is  a threshold  m„(l)  for  Blue  weapons  of  type  i and  for 

D 

each  J = 1,  ...,  N there  is  a threshold  mj^(J)  for  Red  weapons 
of  type  j.  The  battle  terminates  if  B,  (i)  < m„(i)  for  any  i 
or  R^(J)  < m.^(j)  for  any  .1.  That  is,  the  termination  time  T 


This  termination  rule  also  is  subject  to  criticism 
clearly  more  reasonable  than  that  of  (26). 


The  author  begins  with  a heuristic  derivation  of  the  forward 


equation 


for  this  process;  for  a rigorous  derivation  from  a set  of  care 


fully  stated  hypotheses  the  reader  is  referred  to  the  Appendix 
of  [7].  Jennings  proceeds  to  consideration  of  moments  of  first 
and  second  orders,  namely  the  expectations  and  covariances  of 
numbers  of  survivors  given  by 

(28a)  Mg(a;  l,t)  = E“[B^(i)] 

(28b)  Mj^(a;  J,t)  = e“[R^(J)] 

(28c)  Cp(a;i,i’,t)  = e“[.  (B^  (i  )-E“[B^  ( 1)  ] ) (B^  ( 1 • )-E“[B^  ( 1 * ) ] ) ] 

(28d)  Cp(a;J,J',t)  = e“[ (R^ ( J )-e“[ (R^ ( J ) ] ) (R^ (J ' )-E“[R^ (J ' ) ]) ] 
and 

(28e)  Cgp(a;i,J,t)  = e“[(B^(1)-e“[B^(1)]) (R^(J )-E“[R^(j )])]  , 
where  1,1*  =1,  M and  J ,J ' =1,  . . . , N.  In  particular, 

Cg(a;i,l,t)  = Var“(B^(l)) 

and 


Cp(a;j,J,t)  = Var“(R^(J)) 

for  each  1 and  j . 

Jennings  obtains  differential  equations  involving  these 
moments  that  resemble  the  Lanchester  system  of  equations  to 
which  the  stochastic  process  is  analogous.  "Resemble"  in  this 
case  means  that  the  equations  derived  in  [5]  are  qualitatively 
of  the  same  form  as  the  Lanchester  equations,  up  to  an  error 
term  which  in  some  sense  is  small  if  the  time  t is  small. 

For  example,  it  follows  from  Theorem  (21)  that  for  the 
function  f given  by  f(a)  = x^. 
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(29)  it  " 


It  E“[f(B^,R^)] 


= E°‘[Af(B^,R^);.{T>t}] 


r N 
= e“  - I 
L J=i 

( N 

= - E“h 

4=1 


Cj^(J,i)R^(J);{T>t} 


■1 


Cj^(J,l)R^(J) 


] 


r N 1 

+ E“  Cj^(J,i)R^(J);{T<t}J 


N 


= - Cj^(J  ,i)Mj^(aij  ,t) 


+ E 


N 


I Cp(j,l)R. (j);{T<t} 
U = 1 


Similar,  but  more  complicated,  equations  are  valid  for 
variances  and  covariances;  together  -hese  constitute  expression 
(10)  of  Jennings’  paper.  He  gives  no  derivation,  nor  is  there 
explicit  calculation  of  error  terms,  such  as  appears  in  (29); 
we  will  treat  probabilistic  derivation  of  such  equations — 
based  on  Theorem  (21) — and  probabilistic  error  estimation,  in 
a forthcoming  paper  [11]. 

The  author  of  [5]  seems  aware  that  the  error 
A = MB(a,i;t)  - Cb( J ,i )MB(ctiJ  ,t ) j 

is  related  to  the  termination  probability  P°{T<t},  a relation 
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first  observed  in  the  homogeneous  case  by  Snow  [12],  but  seems 
i not  to  know  the  precise  nature  of  the  relation,  which  is  given 

in  (29).  He  proceeds,  nonetheless,  to  attempt  to  compute 
P°‘{T£t}  for  each  t or,  actually,  to  approximate  this  probability. 
From  a computational  standpoint,  approximation  of  these  probabil- 
ities seems  to  us  to  be  of  limited  value.  Such  approximations 
can,  at  best,  warn  one  about  values  of  t for  which  that  approxi- 
mation 

(30)  e“[Bq(1)]  - E“[B^(i)]  = I Cp(J,l)  f E“[R^(J)]du 

.1=1  •'o 

is  grossly  invalid.  While  this  constitutes  useful  Information 
for  very  detailed  combat  models  in  which  the  time  Increment  is 
small,  it  is  precisely  in  such  Instances  that  still  more  care- 
■ ful  error  estimation  is  necessary.  At  the  other  extreme,  for 

highly  aggregated  models  of  large  scale  combat,  it  is  likely 
that  t will  be  sufficiently  large  that  one  incurs  substantial 
errors  by  use  of  (30),  but  no  alternatives  seem  to  exist.  In 
both  cases,  more  accurate  error  estimates  are  required. 

However,  estimates  of  P*^{T£t}  are  better  than  no  error 
estimates  at  all,  so  the  author  deserves  credit  at  least  for 
having  performed  some  preliminary  work.  There  appears  to  be 
* a slight  circularity  to  his  method,  which  consists  in  the 

following  steps: 

(1)  Approximate  Mg(a;i,t)  as  in  (30),  and  perform  similar 
approximations  of  other  moments  defined  in  (28); 

(2)  Assume  that  with  respect  to  P°^,  the  random  vector 
(B^,R. ) is  normally  distributed  with  mean  vector  given  by 

Xf  u 

(28a, b)  and  covariance  matrix  given  by  (28c-e); 


j (3)  Compute  the  mass  of  the  Joint  normal  distribution 

j specified  in  (2)  that  lies  outside  the  set  of  absorbing  states 

1 


(l.e..  In  the  set  of  transient  states  in  which  the  battle  con- 
tinues) and  take  this  as  an  approximation  to  P°‘{T>t}. 

In  truth,  the  author's  approach  is  slightly  different:  he 
fixes  in  advance  an  upper  threshold  value  of  P°‘{T<t}  and  seeks 
the  minimal  t at  which  the  threshold  is  exceeded.  He  effects 
this  computation  by  calculating  the  mass  of  a normal  distribu- 
tion lying  within  an  ellipsoid  on  which  the  Joint  normal  density 
is  constant,  finds  the  time-dependent  constant  density  ellipsoid 
corresponding  to  the  threshold  probability,  and  then  calculates 
the  minimal  time  at  which  the  ellipsoid  intersects  the  termina- 
tion set. 

Virtually  all  these  steps  involve  further  approximations, 
some  of  which  appear  unavoidable  from  a computational  point  of 
view,  at  least  given  the  author's  objective.  The  circularity 
alluded  to  above  is  that  the  assumption  of  a Joint  normal  dis- 
tribution for  (B  ,R, ) is  valid  only  if  p“{T<t}  is  already  very 
small,  and  so  may  not  produce  an  accurate  estimate  of  the  ter- 
mination probability.  Nearly  all  of  the  paper  is  concerned  with 
computations  arising  in  this  approximation  scheme;  the  reader  is 
referred  there  [5]  for  further  details. 
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5.  REVIEW  OF  "STOCHASTIC  'LINEAR  LAW'  BATTLES 

This  paper  [6],  also  by  N.  Jennings,  provides  a treatment 
similar  to  that  accorded  stochastic  "square  law"  battles  in 
the  paper  reviewed  in  Section  2.  The  particular  stochastic 
attrition  process  analyzed  is  Homogeneous  Process  2 of  [10], 
with  proportional  engagement  Initiation  and  single  kills.  The 
infinitesimal  generator  A is  given  by 

A((i,J),  (1,J-1))  = ijcg 

(31)  A((1,J),  (1,J))  = - lj(cg+cp) 

A((1,J),  (i-l,J))  = iJCp  , 

where  c^,  c,,  again  denote  positive  constants,  but  not  with  the 

D K 

same  dimensions  as  the  constants  appearing  in  Sections  2-4. 

The  Jump  function  ^ is  given  by 

X(1,J)  = lJ(Cg+Cj^) 

and  the  transition  matrix  Q of  the  embedded  Markov  chain  is 
given  by 


Q((1,J),  (i,J-l))  = --- Y V 


Q((l,j),  (i-l,J))  = -■ 


In  particular,  the  embedded  Markov  chain  is  a spatially  homo- 
geneous random  walk,  which  is  computationally  tractable. 

Boundary  conditions  in  the  form  of  thresholds  mg  for  Blue 
and  mg  for  Red  are  Imposed. 
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The  author  first  presents  the  forward  equation  for  the 
transition  function  (P^)  of  this  attrition  process.  He  then 
derives  differential  equations  for  moments  of  the  stochastic 
attrition  process  ((B^,R^)).  For  example,  he  presents  the 
equation 


(33) 

Jt 

= - OgE<i"3>[B^R^] 

+ ^[B^R^;{T<t1] 

= - Cj^E^^’J^[B^R^*;{T>t}]  , 


where  T is  the  termination  time  of  the  - engagement . The  author 
provides  only  the  first  equality  In  (33) > we  have  provided  the 
other  two  In  order  to  show  that  this  equation  and  similar  equa- 
tions for  second  moments  can  be  obtained  from  Theorem  (21), 

since  if  (i,j)  is  not  an  absorbing  state  and  f(k,Jl)  = k,  then 

» 

Af (i, j ) = - Cpij  . 

The  author's  equations  appear  correct  (except  for  typo- 
graphical errors)  and  represent  a useful  set  of  facts.  We 
refer  the  reader  to  [6]  for  Jennings'  actual  results. 

The  author  demonstrates  that 
(3M  - c„R,]  - Ogi  - 

for  all  t,  a relation  which  indicates  a plausible  equilibrium 
property  of  this  stochastic  attrition  process.  No  analogous 
property  holds,  however,  for  the  independent  engagement  initi- 
ation process  discussed  in  Sections  2 and  3;  cf.  [11]. 


V-^.' 
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Brief  mention  is  made  of  a computer  program  developed  to 
solve  the  forward  equation  = P^A  numerically;  when  such  a 

solution  is  obtained,  moments  can  be  calculated  at  once.  Of 
course,  this  approach  does  not  lead  to  a general  understand- 
ing of  the  process  in  the  way  that  theoretical  approaches  may. 
Nonetheless,  existence  of  this  program  is  a useful  contribu- 
tion, and.  Indeed,  may  stimulate  and  support  further  theoretical 
developments . 

Finally,  the  author  considers  the  distribution  of  the  termi- 
nation state  (B^,R^).  Let 

q((i,j),(k,il))  = P^^’J^{(B^,R^)  = (k,Jl)}  , 

dependence  of  q on  the  termination  thresholds  nig,mj^  exists  but 
is  suppressed  from  the  notation.  Clearly  q(-,  (k,)!,))  = 0 unless 
k = m„  or  £ = mo,  but  not  both.  At  this  point  the  author's  non- 
probabilistic  approach  leads  him  to  an  unnecessarily  complicated 
derivation  of  this  distribution.  Let 


be  the  probability  that  each  given  casualty  represents  a Blue 

loss.  By  (32)  it  is  immediate  that  for  fixed  initial  conditions 

(1,J),  (Brp,R^)  = (k,mp)  if  and  only  if  the  first  (i-k+( j-mp)-l) 

casualties  represent  i-k  losses  to  Blue  and  (j-mp)-!  losses  to 

st  ^ 

Red  and  the  next  ( ( (l-k)  + (J-mpj) ) ) casualty  is  a loss  to  Red. 

Since  different  casualties  are  Blue  or  Red  with  probabilities 
p and  (1-p),  respectively,  and  are  mutually  Independent  by  vir- 
tue of  the  form  of  the  embedded  Markov  chain,  we  see  at  once 
that 


(35a)  q(  (i,J  ) , (k,nip) ) 


-k  J-nip,-l 

‘"(l-p)  ^ (1-p) 


and  analogously,  that 


(35b)  q( (i, j ) , (mg,£) ) 


(l-m^)  + (j-£)-l\ 


(1-p) 


This  agrees  with  the  result  obtained  by  the  author — his 
equation  (42) — but  provides  more  understanding.  A probabilistic 
derivation  is  illuminating  not  only  in  the  result  obtained  but 
also  in  each  step  and  in  the  overall  pattern  of  reasoning,  in  a 
manner  that  an  analytical  derivation  is  not. 

Let  us  return  now  to  consideration  of  the  terminal  dis- 
tribution. The  distributions  appearing  in  (35)  are  truncated 
negative  binomial  distributions,  further  properties  of  which 
are  discussed,  e.g.,  in  [3]-  Having  derived  the  form  of  the 
terminal  distribution,  the  author  calculates  the  most  probable 
numbers  of  survivors  on  each  side,  given  that  the  other  side  has 
been  forced  to  its  threshold  and  then  considers  a normal  approxi- 
mation to  the  terminal  distribution,  which  appears  to  be  of 
limited  value.  In  any  case,  cf.  [3],  a Poisson  approximation 
may  be  more  appropriate.  We  refer  the  reader  to  Appendix  A of 
[6]  for  details  of  the  results  Jennings  obtains.  Although  these 
results  are  not  discussed  here,  we  do  not  Intend  to  imply  that 
they  are  not  useful  in  the  practical  sense. 


6.  REVIEW  OF  "HOMOGENEOUS  BATTLES  WITH 
GENERAL  ATTRITION  FUNCTIONS" 

This  paper  [15]  by  T.G.  Weale  continues  the  line  of  devel- 
opment represented  by  the  four  papers  reviewed  In  Sections  2-5 
by  extending  some  results  contained  In  the  papers  reviewed  In 
those  sections  to  the  case  of  "general  attrition  functions,”  a 
particular  development  he  pursues  In  [l6,17].  In  our  opinion, 
this  particular  form  of  generalization  has  pitfalls  and.  Indeed, 
can  be  Interpreted  as  an  attitude  which  we  have  long  criticized. 
Except  In  particular  cases  such  as  those  considered  In  [10] 
that  lead,  for  example,  to  the  processes  discussed  In  Sections 
2-5,  there  need  not  exist  a set  of  physical  assumptions  leading 
to  a given  pair  of  attrition  functions.  Mathematical  treatment 
of  processes  not  verlflably  arising  from  physical  assumptions 
about  Individual  weapons  systems  and  their  Interactions  may 
create  the  (possibly)  false  Impression  that  the  resultant  attri- 
tion process  Is  of  value  as  a model  of  combat. 

The  author  treats  a homogeneous  stochastic  attrition  pro- 
cess ((B^,R^))^^q  with  Infinitesimal  generator  A given  by 

A((1,J),(1,J-1))  = «^g(l,J) 

(36)  A((1,J),(1,J))  = - [y,g(l,J)  + V»p(l,J)] 

A((1,J),(1-1,J))  = ^g(l,j)  , 

where  <(>„  and  are  arbitrary  nonnegative  functions  on  the  state 
D n 

space  of  the  attrition  process.  In  addition,  termination  thres- 
holds mg  for  Blue  and  m^  for  Red  are  prescribed  In  the  manner  of 
Sections  2 and  5- 
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That  there  exists  a Markov  process  whose  Infinitesimal 
generator  is  given  by  (36)  is,  of  course,  true.  What  is 
uncertain  is  existence  of  a plausible  or  even  definable  set 
of  assumptions  concerning  physical  behavior  of  combatants — 
both  individually  and  interactively — that  leads  to  an  attrition 
process  with  the  generator  A given  by  (36).  For  certain  cases, 
the  author  of  this  review  has  shown  the  existence  of  such  sets 
of  assumptions;  cf.  [7jl0].  These  assumptions  concern  the  quali- 
tative and  quantitative  probabilistic  nature  of  engagement 
initiation  by  combatants  and  are  stated  in  a form  in  which  a 
potential  appller  of  stochastic  attrition  processes  (to  com- 
puterized combat  simulations,  for  example)  can  readily  verify 
their  plausibility,  or  at  least  choose  one  process  among  several 
alternatives. 


Even  from  a mathematical  standpoint,  the  arbitrary  generator 
approach  is  subject  to  criticism.  Except  for  the  restriction 
that  sample  paths  be  componentwise  nonincreasing  and  decrease 
only  by  Jumps  of  size  one  in  one  component,  the  infinitesimal 
generator  A given  by  (36)  is  perfectly  general.  Specific  com- 
putations are,  therefore,  likely  to  be  impossible  to  perform, 
as  the  main  body  of  [15]  confirms.  Theoretical  results  of 
sufficient  specificity  to  be  of  interest  are  likewise  difficult 
to  obtain.  The  role  of  (restrictive)  assumptions  in  mathematics 
is  to  sufficiently  limit  the  class  of  objects  under  study  that 
nontrivial  statements  become  possible.  Furthermore,  considera- 
tion of  problems  that  are  too  general  denies  one  use  of  both 
intuition  and  methods  of  analysis  that  exploit  the  special 
structure  of  specific  problems. 

Finally,  one  does  not  deal  in  practice  with  an  "arbitrary" 
generator  A but  chooses  some  specific  form.  A set  of  physical 
assumptions  leading  to  a given  form  of  generator  contains  with- 
in Itself  the  appropriate  suggestions  for  generalization.  If  a 
certain  assumption  is  believed  to  be  implausible,  one  can  modify 
it  to  be  acceptable  in  physical  terms  and  then  derive — possibly 
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not  without  difficulty — the  generator  resulting  from  the  new 
set  of  assumptions.  On  the  other  hand,  a methodology  based  on 
direct  and  arbitrary  choice  of  the  generator  Is  Inherently  self- 
limiting  in  that  it  admits  no  such  potential  for  generalization. 

For  all  these  reasons,  we  believe  that  the  research  re- 
ported by  Weale  in  this  paper  has,  in  its  present  form,  little 
significant  practical  or  mathematical  implication.  Therefore, 
the  following  description  of  the  contents  of  the  paper  is  quite 
brief. 

Weale  takes  note  of  the  forward  equation 
(37)  = P^A 

with  suitable  boundary  conditions,  for  the  transition  function 
(P^)  of  the  attrition  process.  No  derivation  is  required,  for 
Markov  process  theory  ensures  the  validity  of  the  forward  equa- 
tion. For  reasons  that  are  obscure  to  us,  the  author  of  [15] 
thereafter  considers  higher  derivatives  of  the  transition  func- 
tion and  obtains  various  relations  involving  them.  It  is 
evident  from  (37)  that 


P"  - P'A  - P^A- 


and  that,  more  generally,  for  each  n 
(38)  P^"^  = 

for  all  t and  k = 0,  ...,  n,  where  is  the  k^^  derivative 

of  the  transition  function  and,  by  convention  P^^^  = P^ • The 
utility  of  (38)  in  the  context  of  stochastic  attrition  processes 
is  not  apparent  to  us,  as  (37)  is  quite  sufficient  to  specify 
(P. ) uniquely.  For  general  attrition  functions  further  Infor- 
matlon  given  by  (38)  is  not  sufficiently  specific  to  be  of 
value.  In  special  cases,  of  course,  a useful  Taylor  expansion 
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of  might  be  obtainable,  but  the  expansion  follows  more  easily 
from  the  well-known  property  that  = exp  (tA). 

The  author  Includes  some  remarks  concerning  the  qualitative 
form  of  the  attrition  functions  For  example,  both  should 

be  increasing  in  each  variable  separately  and  (possibly)  strictly 
increasing  in  the  variable  representing  the  opposition.  For 
example  <Ppj(i,J),  which  corresponds  to  Red  kills  of  Blue  weapons 
should  be  nondecreasing  in  1 (which  represents  the  Blue  side)  and 
strictly  increasing  in  J.  The  author  proposes  that  J -►  «^>o(i,j) 
be  a function  with  an  S-shaped  graph,  a property  not  possessed  by 
any  of  the  attrition  functions  so  far  derived  from  physical 
assumptions  (which  are  linear  in  j ) , but  certainly  plausible. 

A discussion  is  given  of  the  general  form  of  the  solution 
to  (37);  cf.  Section  2 for  a similar  treatment  of  a specific 
case.  Even  in  that  specific  case  no  concrete  results  are  ob- 
tained; in  this  general  case  nothing  of  Interest  is  presented. 

Weale  mentions  in  the  main  text,  and  presents  in  Appendices, 
computer  programs  for  numerical  integration  of  the  forward  equa- 
tion (37)  together  with  some  results  obtained  therefrom.  The 
attrition  functions  used  in  the  sample  program  are  given  by 

- KV'b'” 

and 

V>j^(i,j)  = j(Cp+C^i) 

where  c_,  c',  c„,  c*  are  constants.  We  emphasize  that  to  our 

D D R n 

knowledge  no  set  of  physical  assumptions  leading  to  these  attri- 
tion functions  exists.  That  such  assumptions  do  exist  is. 

Indeed,  quite  possible;  verification  of  such  matters  is  the 
problem  of  Interest  and  importance. 
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Finally,  the  author  treats  the  terminal  distribution  by 
the  same  Tauberlan  methodology  used  in  the  paper  reviewed  in 
Section  2.  No  specific  analytical  results  are  obtained, 

although  many  numerical  explorations  are  possible.  i 
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7.  REVIEW  OF  "THE  DISTRIBUTION  OF  THE 
DURATION  OF  BATTLE" 

In  this  paper  [l6]  T.G.  Weale  continues  his  analysis  of 
the  homogeneous  attrition  process  ((B^,R^))^^q  introduced  in 
[15]  (cf.  Section  6 for  a review  thereof):  namely,  the  Markov 
attrition  process  with  infinitesimal  generator  A given  by 

A((i,J),(i,j-l))  = -TgCi.j) 

(39)  A((i,j),(i,J))  = - [c^.3(i,j)  + 'Tj^(i,J)] 

A((i,j),(i-l,J))  = ^fp(i,j)  , 

where  cp  cp  are  arbitrary,  but  fixed,  nonnegative  functions. 

D n 

As  we  have  observed  in  Section  6,  this  general  approach  has 

difficulties:  for  most  functions  cp_,cp  there  is  no  known 

D n 

set  of  underlying  physical  assumptions  that  leads  in  the  manner 
of  [7,10]  to  an  attrition  process  with  the  generator  A given  by 
(39).  Moreover  the  general  approach  represented  by  [15,16,17] 
appears  unlikely  to  yield  results  of  sufficient  specificity  to 
be  really  useful. 

Nonetheless  the  particular  problem  studied  in  [I6]  is  of 
some  interest — especially  in  physically  justifiable  special 
cases  such  as  the  processes  discussed  in  the  DOAE  papers  [131, 
[14],  and  [6]  that  are  reviewed  in  Sections  2,  3>  5,  respectively. 
That  problem  is  the  following:  let  mg  and  m^  be  termination 
levels  (possibly  but  not  necessarily  zero)  for  the  Blue  and  Red 
sides,  respectively,  so  that  the  combat  ceases  when  B^  = mg  nr 
R.  = m„.  Let  T be  the  duration  of  the  battle,  given  by 

Z A 

T = lnf{t:  B.  = m„  or  R.  = m„}  . 

U b U A 
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One  then  wishes  to  compute,  characterize,  or  approximate,  for 
each  Initial  state  (i,j),  the  distribution  of  T under  the  proba- 
bility measure  l.e.,  to  study  the  function 

(40)  P(i,j;t)  = < t}  . 

As  observed  in  Section  4 (cf.  page  21)  such  probabilities  are 
of  interest  in  the  context  of  the  moment  equations  derived  in 
Theorem  (21)  in  Section. 3 and  are  also,  of  course,  of  Intrinsic 
interest  as  properties  of  the  combat. 

Weale's  approach  to  the  problem  is  the  following:  if 
denotes  the  state  of  the  process  at  the  time  of  termination 
and  if  1 > m„,  j > m_,  then 

D K 

€ D}  = 1 , 

where 


D = {(k,mp):  m^  < k}  u {(m^.A):  m^^  < £}  . 

Moreover,  since  each  state  in  D is  by  definition  an  absorbing 
state,  it  follows  that 

(41)  F(l,J;t)  = C D} 

i 3 

= V P,  ((i,j),(k,m„))  + 2^  P.((i J),(m„,Jl)) 

k=m^+l  ^ ^ S,=mp+1  ^ ^ 

The  expression  given  in  (41)  is  valid  in  general,  whether  the 
generator  A arises  from  explicit  physical  assumptions  or  not. 
Therefore,  if  the  transition  function  ( P^ ) of  the  process  were 
known,  then  the  functions  P(i,j;*)  would  be  completely  deter- 
mined. Unfortunately,  however,  this  is  not  so  even  for  the 
special  cases  discussed  in  previous  sections. 
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Weale  partitions  the  termination  set  D into  three  subsets 
consisting  of 


= { (k,mpj) : k > nig}  , 

D2  = {(mg,)l)  : I > nig} 

and 

= D - (D^  U D^)  , 

where  m_  > m„  and  m_  > m„  are  also  prescribed  in  advance.  States 

D — D n — K 

in  are  in  Weale 's  terminology  "Blue  victory  states,"  those 
in  are  "Red  victory  states,"  and  those  in  are  "draw  states. 
The  intention  and  Interpretations  are  apparent.  Clearly 

p(l,J){B^  C D^}  = P„((i,j),(k,mg))  , 

k-mB 


(42) 


^{B^  C D2} 


i P<„((i,J),(mg,fc))  , 


P^^’^^{B„  C D„}  = I P^((l,J),(k,m  ) 

^ ^ k=m„+l  “ 

D 

niR-l 

+ I P„(  (ij  ) , (niD.i) ) 
il=m„+l  “ 

n 


= 1 _ p^^’J C U D2}  , 

where 


P 


00 


11m  P, 
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mu: 


and  evidently  exists  provided  that  <Pn(i,J)  + <P„(1,J)  > 0 when- 
ever  1 # nig  and  j (cf.  pp.  8-9  above  for  some  remarks  rele- 

vant to  numerical  computation  of  P^) • 

Weale  then  proceeds  to  consider  the  distribution  of  T 
conditioned  on  the  termination  state;  that  is,  he  considers 
the  conditional  distribution  functions 

F^(i,j;t)  = < t)B^  C D^} 

F^CiJ;!)  = < tlB^  C D^} 

F^(i,j;t)  = < t|B^  C D^}  . 

Formally,  these  functions  can  be  computed  directly  from  (4l) 
and  (42)  using  only  elementary  probability;  for  example 

p(i,j)rTj  _ n 1 


UB.  c Dt  } 
1 p(i,J c D, } 


P^(  (1,J  ) , (k,m  ) ) 

k^mg 

i 

k=mB 

Analogous  expressions,  which  need  not  be  included  here,  exist 
for  p2(l,J;t)  and  F2(l,j;t).  From  (43)  it  is  clear  that  the 
Fq(i,J;t)  can  be  computed  in  closed  form  if  the  transition 
function  (P^)  is  known,  but  seem  inaccessible  otherwise. 

To  continue  our  review  of  [16],  Weale  next  discusses  the 
density  functions  of  the  distributions  P(l,j;t).  From  (4l)  and 
the  forward  equation  Weale  derives  the  relation 


(4H)  ^F(l,j;t) 


1 J 


B 


‘R 


= y P.  A((i,J),(k,inj,))  + V P.A((iJ),(m„,£)) 
k=f^+l^  ^ ^ 


= ^ P^((l,J),(k,mp-rl))q''g(k,mj^+l) 


k=mg+l 


O 

+ P^((i,j),(nig+l,«,))crg(mg+l,Jl)  . 


il=mg+l 


Analogous  expressions  are  also  given  for  the  densities  of  the 
conditional  distributions  F^,  F^,  Like  other  relations 

derived  earlier  in  the  paper  these  formulas  give  closed  form 
results  if  the  transition  function  (P^)  is  known  in  closed 
form  and  not  much  information  in  other  cases. 

One  exception  to  the  latter  assertion  is  when  the  forward 
equation  is  Integrated  numerically  for  given  initial  conditions 
and  attrition  functions;  this  is  a principal  objective  of  the 
computer  programs  that  comprise  the  larger  part  of  not  only  [l6] 
but  also  several  of  the  papers  reviewed  in  preceding  sections  of 
this  paper.  In  this  situation  the  formulas  (^1),  (42),  (43), 
and  (44)  are  all  applicable  and  yield  good  numerical  approxi- 
mations to  the  probabilities  in  question.  The  work  and  results 
concerning  numerical  computations  are,  as  we  discuss  further  in 
Section  9,  a significant  contribution  of  at  least  some  of  these 
DOAE  papers.  This  method  of  development  for  all  its  usefulness, 
however,  does  not  yield  closed-form  analytical  expressions  by 
means  of  which  one  can  fully  understand  and  describe  the  attri- 
tion processes  under  study.  Neither  does  it  produce  qualitative 
insights  that  can  be  extrapolated  beyond  the  (necessarily 
limited)  numbers  of  numerical  inputs  that  are  actually  treated; 
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of  course  many  useful  qualitative  insights  result  nevertheless. 

Finally,  Weale  deals  with  expectations,  medians,  and  modes 
of  the  termination  time  distributions.  Since,  for  example,  for 
each  (i,j) 


(^5) 


(l-F(i,J;t))dt 
(1-F(1,J ;t))dt 


u 


u 

F(i,J;t)dt 


for  large  values  of  u,  one  can  use  the  first  equality  in  (45) 
to  compute  if  F(1,J;*)  is  known  in  closed  form  or  can 

use  the  approximation  contained  in  (45)  if  (as  is  the  case  in 
[l6])a  numerical  approxiraation  to  F(i,J;-)  is  available.  Since 


u 


< u. 


in  order  that  the  approximation  in  (45)  be  even  reasonably 
accurate  one  must  have 

u > ^[T]  . 


This  point  is  not  mentioned  in  [l6];  it  is  easy,  however,  to 
estimate  E^^’'^^[T].  If,  as  the  physics  of  combat  almost  cer- 
tainly require,  the  attrition  functions  '!■  are  both  nonin- 

D n 

creasing  in  each  argument,  then  evidently 

(46)  < (l-mg+J-mj^-l)[o.g(mg+l,mp+l)  + rp(mg+l,mp+l)]"^ 


The  reasoning  underlying  (46)  is  the  following:  at  most  (i-m- 

D 

+ J-nip  - 1)  casualties  can  occur  before  termination  and  each 
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interval  between  casualties  Is  exponentially  distributed  with 
expectation  not  exceeding  [cp 

D 

One  may  then  use  (46)  to  determine  a sufficiently  large  value 
of  u for  use  In  (45).  In  the  computer  program  appended  to 
[16]  the  integration  In  (45)  is  performed  using  Simpson's  rule, 
which  seems  eminently  reasonable. 

Computation  of  the  median  and  mode  of  the  distribution  of 
T is  also  considered  in  [l6].  Also,  the  author  provides  simi- 
lar treatments  of  the  conditional  distributions  F^,  P^,  P^J 
further  details  are  not  necessary  here. 

As  we  mentioned  before,  the  distribution  of  the  terminal 
time  T is  of  some  interest,  particularly  in  the  context  of  the 
differential  equations  discussed  in  Section  3;  see  also  [11] 
where  we  consider  the  problem  in  some  detail.  One  approach 
that  may  be  superior  to  that  of  Weale  not  only  for  derivation 
of  analytical  results  but  also  for  certain  computational  appli- 
cations is  the  following  recursive  method. 


(mg+l,mp+l)  + cpp(mg+l,mj^+l)] 


(47)  PROPOSITION.  For  each  nonabsorbing  state  (i,j) 


(48) 


P(i,J;t) 


/.t  -C<Pn(i,J  ) + 4'R(l,j  )]u 

= ‘VgCl,])  j F(i,J-l;t-u)e  ^ ^ du 

0 


F(i-l,j ;t-u)e 


-[TB(i.J)+'^R(i  J)]u 


du 


PROOF.  Let  be  the  first  change  of  state  of  the  attrition 
process  ((B^,R^))  and  the  state  entered  at  time  T^.  By 
Theorem  (8.3.3)  of  [1], 
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cpg(l,J  )+cfg(i,J  ) 

-[(;p  (i,J)+4'  (l,J)]u 

[^B(i,J)+crR(l,J)]e  ® ^ du 


t )+cpp(i,J  )]u 

= e ® ^ du 


and,  in  the  same  way. 


p(l,J){T^<t;X^=(i-l,j)}  = 'Pp(i,j)  J 


t -[cpg(i,J)+Cf'  (i,J)]u 
e du  . 


Hence, 

F(l,J;t)  = 

= E^^’‘^^[P  ^{T<^t-Tj^}] 

= E^^’*^^[F(X^;t-T^)] 

- CD  ri  n -[q'  (i,j)+i^;g(l,J)]u 

- VgdJ)  / F(i,j.i;t-u)e  ® ^ du 

-['1'  (l,j  )+VR(i,J  )]u 

++p(l,J)  f F(i-l,J;t-u)e  ^ " du  , 

*^0 

where  the  third  equality  Is  by  the  strong  Markov  property; 
cf.  [1]  or  [2].  0 


J 


) 


To  illustrate,  we  consider  the  process  with 


¥>g(l,j)  = Cgi 


= Cj^J  , 


namely  the  homogeneous  square  law  process  discussed  in  Sections 
2 and  3 above.  In  this  and  other  applications  it  is  more  con- 
venient to  use  (48)  in  the  equivalent  form 

(49)  p^^><JMT>t}  = vJgdJ)  f ■^^{T>t-u}e  ® ^ c 

•'o 

-t  , .X  -[^  (i,J)  + v>p(i,j)]u 

+ v’p(l,J)  f nT>t-u}e  ® ^ dv 

) + ^p(l,J  )]t 

^6  • 

For  < 2 we  then  have  the  following  exact  results: 

p(0>0)(ip>t}  = P^^»°^{T>t} 

= P^°’^^{T>t} 

= P^^’°^{T>t} 

= p(°»2){^>t} 

= 0 , 


while 


Also, 


-(Cb+Cr)! 


p(ljl){,p>t  } = e 


/p  -I  \ Cp  “(Cp+Cp)t  Cp  — (2cQ'^Cp)t 

P^^’^^{T>t}  = ^ e ^ ^ + (1— ^)e  ^ ^ 


i 

i 


8.  REVIEW  OF  "MOMENTS  OF  THE  DISTRIBUTION  OF  STATES 
FOR  A BATTLE  WITH  GENERAL  ATTRITION  FUNCTIONS" 

This  paper  [17]  by  T.G.  Weale  and  E.  Peryer  continues  the 
work  of  the  first  author  on  homogeneous  battles  with  general 
attrition  functions  that  is  reported  in  [15]  and  [16],  and  are 
reviewed  in  Sections  6 and  7 above,  respectively.  The  general 
comments  made  in  those  sections  (in  Section  6 in  particular) 
concerning  possible  lack  of  an  underlying  family  of  physical 
assumptions  remain  relevant.  In  [17]  Weale  and  Peryer  deal 
with  results  analogous  to  those  obtained  in  [5]  by  Jennings 
for  the  heterogeneous  square  law  attrition  process;  namely, 
differential  equations  for  expectations  of  functionals  of  the 
attrition  process.  Many  of  the  comments  and  analyses  presented 
in  Section  4 of  this  review,  in  which  [5]  is  reviewed,  will 
also  be  germane  to  the  discussion  of  [17]. 

For  the  sake  of  completeness  we  once  again  observe  that 
the  attrition  process  treated  in  [17]  is  a Markov  process 
( (B^  ,R^ ) with  infinitesimal  generator  A given  by 

A((i,J);(l,J-l))  = 

A((l,J);(i-l,j))  = , 

where  and  <f>„  are  nonnegative  but  otherwise  arbitrary  functions 
defined  on  the  state  space  E = N x N of  the  attrition  process. 

The  authors  first  take  note  of  the  forward  equation  for  the 
transition  function  of  the  process,  which  in  open  form  is  given 
for  nonabsorbing  states  (k,Jl)  by 
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(50)  P^((l,J),(k,Jl))  = P^A((l,J),(k,]l)) 

= Pt((l,J),(k,il+l))^g(k,il+l) 

+ Pt((l,J),(k+l,ll))^p(k+l,Jl)  , 


for  absorbing  states  (k,mpj),  where  Is  the  Red  termination 


level,  by 


(51a)  P(((i,J),(k,mj^))  = P^(  (1,J  ) , (k.m^+l)  )v»g(k,mj^+l)  , 


and  for  absorbing  states  (mg,jl),  where  is  the  Blue  termina- 


B 


tlon  level,  by 


(51b)  P^((i,J),(mg,)l))  = P^((l,J),(mg+l,£))v»j^(mg+l,Jl). 


It  then  follows  that  if  f is  a function  on  the  state  space 
E of  the  attrition  process, 

i .1 

(52: 


I 


it  “ it  23  S f(k,£)p  ((i,j),(k,ji))! 

k=m„  )l=m„  * 


B R 


= X]  f(k,£)P'((l,J),(k,)l)) 


k=mg 


2]  f(k,)l)[v^(k,£+l)P^((l,J),(k,Jl+l)) 


k=m_+l  £=mD+l 
B n 


+ ^g(k+l,£)P^((i,j),(k+l,£))] 


+ 53  f(k,mg)^g(K,mg+l)P^((i,J),(k,mg+l)) 

k=mg+l 

J 


+ V f(mg,«)^p(ing+l,£)P^((l,1),(mg+l,ll)) 


£=f^l 


The  expression  (52)  is  essentially  identical  to  equation  (8)  in 
[17]  which  is  the  main  theoretical  result  therein.  Observe  that 
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the  right-hand  side  of  (52)  Involves  no  probabilities  of  the 
forms  P^(  (1,J  ) , (k,mp^) ) or  ( ( i , j ) , (mg, H) ) . Indeed,  elementary 

calculations  verify  that  (52)  Is  equivalent  to  the  result  that 
Theorem  (21) — as  suitably  extended  in  [11]  to  apply  to  general 
Markov  attrition  processes — would  yield  in  this  situation. 

The  authors  then  proceed  to  discuss  the  solution  of  (52) 
for  functions  f of  the  form 

f(k,«,)  = (k-l)(k-l+l)  •••  (k-l+r)()l-J ) •••  (Jl-j+s)  , 

in  which  case  the  expectations  involved  are  factorial  moments. 
They  further  assume  that  and  are  polynomials  in  their 

D n 

arguments,  but  even  so  do  not  obtain  specific  results.  The 

difficulty,  the  reviewer  believes,  is  that  the  authors*  point 

of  view  and  method  of  proceeding  are  analytic  rather  than 

physical  and  probabilistic.  From  an  analytic  standpoint  there 

seems  to  be  hope  of  solving  (52)  only  if  the  and  are  poly- 

D n 

nomlals,  although  useful  approximations  involving  polynomials 
can  certainly  be  made.  From  the  physical,  probabilistic  stand- 
point one  should  attempt  to  solve  (52)  only  for  attrition 
functions  arising  from  well-defined  physical  hypotheses,  in 
which  case  the  particular  probabilistic  structure  of  the  process 
at  hand  may  aid  in  obtaining  a solution. 

Appendices  to  the  paper  describe  computer  programs  designed 
to  approximate  the  solution  of  (52).  There  is  no  doubt  that  one 
can  solve  (52)  numerically  for  essentially  any  attrition  func- 
tions that  can  be  programmed  into  a computer,  although  for 
irregular  functions  the  approximation  may  not  be  good.  For 
attrition  functions  verifiably  arising  from  physically  definable 
and  plausible  assumptions,  these  computer  programs  constitute  an 
important  analytical  and  descriptive  tool. 


9.  CONCLUSIONS 

The  first  four  papers  presented  here  represent  varying  but 
substantive  contributions  to  the  theory  of  stochastic  attrition 
processes.  To  the  taste  of  the  reviewer  these  contributions 
are  lessened  and  obscured  by  the  authors*  excessive  reliance  on 
analytic  approaches  since,  as  we  have  several  places  demonstrated, 
direct  probabilistic  approaches  yield  not  only  more  Illuminating 
arguments  but  also,  at  least  In  some  cases,  more  complete  or 
specific  results.  The  computer  programs  associated  with  these 
papers  constitute  a contribution  whose  current  value.  In  view  of 
our  lack  of  ability  to  deal  with  the  processes  on  a reasonable 
closed-form  basis.  Is  certainly  understated  by  the  scant  atten- 
tion devoted  to  them  In  this  review.  In  a few  places  lengthy 
and  unenllghtenlng  computations  Interrupt  the  development  of 
worthwhile  Ideas,  but  In  general  the  papers  are  of  high  techni- 
cal and  expository  quality,  for  which  their  respective  authors 
are  to  be  commended. 

On  the  other  hand  we  believe  that  the  last  three  papers 
dealing  with  battles  with  general  attrition  functions  make  little 
contribution  In  the  mathematical  sense.  Since  our  criticisms  are 
rather  strongly  worded,  let  us  once  more  attempt  to  be  specific 
about  our  grounds  for  criticism.  The  first  of  these  Is  philo- 
sophical but  has  Important  practical  Implications;  we  strongly 
believe  that  one  should  not  deal  with  attrition  processes  that 
cannot  be  Justified  In  terms  of  physical  assumptions. 

But  our  criticisms  are  also  on  mathematical  grounds:  the 
great  degree  of  generality  Involved  prohibits  the  authors' 
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obtaining  results  that  are  sufficiently  specific  to  be  of  mathe- 

matical  interest  even  without  regard  to  possible  physical  Interest  | 

or  applicability.  | | 

Of  much  greater  Importance  and  in  much  greater  need,  ir  our 
opinion,  is  work  aimed  at  developing  physical  assumptions  that 
imply  forms  of  attrition  functions  other  than  those  previously  ! j 

(in  [10],  e.g.)  justified.  Better  yet,  one  should  strive  first  \\ 

to  develop  plausible  sets  of  physical  assumptions  and  then  to  \ ! 

derive  from  these  attrition  processes  that  are  tractable  in  terms  ’ | ; 
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